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(54) Soft gain update method & apparatus 

(57) This invention is concerned with updating (or 
changing) volume (or gain) in a soft manner for voice/ 
speech and 1 0% audio applications. It is found that, from 
mute to full-scale gain update, filtered gain update pro- 
cedure (in line or domain) provides subjectively superior 
performance over instantaneous or zero-cross gain up- 
date or step-based update in the decibel domain. Spe- 



cifically, it is found that about 90 milliseconds (ms) of 
update time for first-order (linear) update and about 45 
milliseconds of update time for second-order (quadratic) 
update provide the least audible distortion with the 
quadratic update performing significantly better than the 
linear first-order update. Both filtered updates are car- 
ried out in the linear domain, not in the decibel domain. 
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Description 

Field of the invention 

5 [0001] The present invention relates to updating of gain in audio amplification, particularly to updating of gain in 
complex integrated circuits known as digital signal processors (or "DSP"). 

Background of the Invention 

w [0002] Traditionally, gain or volume in an audio context is updated instantaneously when an update command is 
received by a gain processor. However, when the gain is changed from its present value of G old to a new value, G new , 
the audio signal exhibits an abrupt jump (or "discontinuity") at the transition instance. A gam update can also be delayed 
until the signal to which the volume change will be applied is crossing the zero value. This is "zero-cross" gain update. 
Zero-cross gains usually operate in a time-out mode which forces a gain update if a zero-cross of the signal does not 

'5 occur within a prescribed time, about 1 00 ms. For a constant or dc signal, both instantaneous and zero-cross updates 
yield identical results. 

[0003] Figure 3 depicts this time-domain phenomenon due to a gain signal. g(n), being multiplied by a constant data 
signal, x(n) = c. The gain applied signal, y(n) = x(n)g(n), results. Discontinuity at time 0 is clearly exhibited. 
[0004] With a dc signal, an instantaneous gain-update causes a sharp transient whose energy is distributed over 
20 the entire audible band (0 < f < fs/2), causing unpleasant audible effects. Figure 4 depicts this phenomenon in the 
frequency-domain due to a gain signal, G(/), where 0 < f< fs/2, being multiplied by a constant data signal x(0 = c5(0- 
the discrete-time fourier transform of x(n). The gain applied signal: 

vm - G(frx(f) - G ( f )* cb (0 - £9Q 
25 [) ~ 2n ~ 2n ~ 2k 

results. 

[0005] However, most signals are non-constant in that they are time-varying. Prior hardware-based solutions have 
employed zero-cross gain update methodology, which provides slightly better audible performance than the instanta- 
30 neous gain update method for non-constant signals. The application of new gain to non-constant signals during zero 
crossing of x(n) helps reduce the audible artifacts to a small degree. 

[0006] Figure 5 depicts the results of instantaneous and zero-cross gain updating. A gain update command 502 is 
received simultaneously for both instantaneous 504 and zero-cross gain update 506. For an instantaneous gain update, 
the gain is updated at time 508 when the command is received. This immediate update results in a signal spike 510 
is causing audible audio artifacts. For a zero-cross update, the gain is updated after the command is received, once the 
signal is near the zero value 512. However, the signal is stilt updated instantaneously at the zero-cross, again resulting 
in audible audio artifacts at time 514. 

[0007] A "soft mute" is another current gain update implementation. In soft mute gain update, the gain level is grad- 
ually reduced to zero (-<*> in the decibel domain) over a time frame that is easily perceptible by the listener (e.g. several 
40 seconds). This implementation serves to reduce the audible artifacts created by instantaneous or zero-cross updates. 
However, the gain update time is detectable by the listener. 

Summary of the Invention 

45 [0008] The present application describes an automatic soft gain update architecture in which the gain is automatically 
swept from value G otd to value G new in the linear domain over a period of time which is imperceptible to the listener. 
Such a gain sweep reduces the unpleasant audible effects caused by discontinuity during instantaneous or zero-cross 
updates which take place in the decibel domain and eliminates the perception of the length of the update by the user. 
Two specific soft gam update implementations, namely first-order (or "linear"; update and second-order (or "quadratic") 

50 update will be described. 

[0009] The first-order gain function is continuous. However its derivative is discontinuous at 0 and N v Linear gain 
update can be implemented in either a constant-step or constant-time approach. In the constant-step approach, the 
gam is automatically swept linearly over a number of steps determined by G oid and G new . The number of steps deter- 
mines the time taken for the gam sweep to complete. In the constant-time approach, the gain is automatically swept 

55 linearly from G okJ to G new . Over a specified period of time, regardless of the.number of steps between G otd and G new . 
Sweep time for linear gain updates is , where T, = Ts, with N, indicating the number of steps required for the gain 
sweep and Ts being the sampling period. 

[0010] The second-order gain function and its first derivative are both continuous. However, its second derivative is 
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discontinuous. In the quadratic constant-time embodiment, the gain is automatically swept according to a ramping 
function. The duration of time of the gain sweep from G o!d to G new is of significant importance as well. Sweep time for 
a quadratic gain update is T 2 , where T 2 = N 2 Ts, with N 2 indicating the number of steps required for the gain sweep 
and Ts being the sampling period. 

5 [0011] One advantage of the present disclosure is that polynomial gain update procedures allow significant reduction 
in audible artifacts. Another advantage of the present disclosure is that the polynomial gain update can be performed 
in a linear (first-order), quadratic (second-order), or higher-ordered fashion. Another advantage is that from mute to 
full-scale gam update, a filtered gain update procedure in the linear domain, provides subjectively superior performance 
over instantaneous or zero-cross gain update or step-based update in the decibel domain. Another advantage is the 

10 guarantee of perfect gain update in that there are no gam errors or approximations. That is, the gain value G new is 
achieved exactly as opposed to an analog representation which may deviate from G new by 2% or more. 

Bfief Description of the Drawings 

'5 [0012] The present invention will now be further described, by way of example, with reference to the accompanying 
drawings, which show important exemplary embodiments of the present invention wherein: 

Figure 1 represents the quadratic constant-time filter soft gain update: 

Figure 2 depicts a block diagram of the signal processing of an audio system which utilizes the soft gain update; 
20 Figure 3 depicts the time-domain phenomenon due to a gain signal being multiplied by a constant data signal; 

Figure 4 depicts the frequency -domain phenomenon due to a gain signal being multiplied by a constant data signal: 

Figure 5 depicts the results of instantaneous and zero-cross gain updating; 

Figure 6 depicts linear soft gain update response in the decibel domain; 

Figure 7 depicts linear soft gain update response in the linear domain: 
25 Figure 8 depicts the response of a constant-step soft gain update; 

Figure 9 depicts the response of a gain update from 1/2 full scale to full-scale; 

Figure 10 depicts in flowchart form the algorithm for controlling a linear constant-step gain update process; 
Figure 11 illustrates a gain processor which implements a constant-step linear sweep gain update; 
Figure 12A depicts the constant nature of in a constant time gain update which sweeps from mute to maxgain; 
30 Figure 12B depicts the constant nature of N., in a constant time gain update which sweeps from maxgain/2 to 

maxgain, 

Figure 13 represents the linear constant-time filter output desired; 

Figure 14 schematically depicts the implementation of a linear constant-time difference equation: 
Figure 15 illustrates a gain processor which implements a quadratic constant-time sweep gain update: and 
35 Figure 16 depicts a microcode cycle-based DSP implementation of a table algorithm with a single adder. 

Detailed Description of the Preferred Embodiments 

[0013] The numerous innovative teachings of the present application will be described with particular reference to 
-to the presently preferred embodiment. However, it should be understood that this class of embodiments provides only 
a few examples of the many advantageous uses of the innovative teachings herein. In general, statements made in 
the specification of the present application do not necessarily delimit any of the various claimed inventions. Moreover, 
some statements may apply to some inventive features but not to others. 

[0014] The present application provides a system of gain update in which the gain is automatically swept from G 0 , d 
•ts to G new over a period of time in order to produce a gain update with is audibly smooth. In various embodiments, the 
period of time over which the gain is swept is derived in two ways. The first method described below is the constant 
step method. In the constant step method, the range of gain between mute and full scale is divided into steps. The 
time required to update gain is the same for each step. Therefore, gain updates with large differences between G oW 
and G new take longer to complete than those with smaller differences. The response of a linear constant-step soft gain 
50 update is depicted in Figure 8. Linear constant-step soft gain update achieves a reduction in audible artifacts at gain 
update. 

[0015] The second method described below is the constant-time method. In the constant time method, the time 
required to update gain from G old and G new is constant for each update regardless of the difference between G old and 
G new . The response of a linear constant-time soft gain update is depicted in Figures 1 2A and 1 2B. Like linear constant- 
's step updating, a linear constant-time soft gain update achieves a reduction in audible artifacts. However, the quadratic 
constant-time soft gain update, as depicted in Figure 1 , results in a significant further reduction of audible artifacts due 
to its first-order continuity at G new . 

[0016] When sweeping from mute to largest gain, sweep times of about 90ms for linear gain updates and about 
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45ms for quadratic gain updates are found to have the most pleasant listening experience. In a 1 6-bit audio application, 
for instance, the greatest gain update sweep is from mute, y(n) = 0, to full-scale, ly(n)l = 2 1S - 1 = 32667. Values of Nj 
and Tj, i = 1 or 2, for sampling rates of 44 1 kHz and 48 kHz, respectively, are shown below: 



Fs 


Ni 


Ti 


N 2 


T 2 


Ts (sampling period) 


44.1 kHz 


4096 


92.88 ms 


2048 


46.44 ms 


22.68 U-S 


48.0 kHz 


4096 


85.33 ms 


2048 


42.67 ms 


20.83 u,s 



io [0017] Linear sweeping of the gain from one value to another has been implemented in software and in hardware 
However, such sweep functions typically take place over a non-optimum period of time or in the decibel domain Such 
sweep functions cause audible distortion. About 90 milliseconds (ms) of update time for first-order (linear) updates and 
about 45 milliseconds of update time for second-order (quadratic) updates provides the least audible distortion due to 
sweep time, with the quadratic update performing significantly better than the linear first-order update. 

is [001 8] 1 1 is also important, in the presently preferred embodiment, that gain sweeping take place in the linear domain, 
not m the decibel domain, where g dB (n) = 20log 10 (g (n) ), g (n) > 0. If gain is swept in constant dB steps, e.g. from -90dB 
to OdB in steps of 1 dB, a decibel response like that depicted in Figure 6 results. 

[0019] The linear-domain response equivalent of Figure 6 is depicted in Figure 7. Applying a constant step-update 
t° 9dB( n ) translates to exponential gain update in the linear-domain gain, resulting in significant audible artifacts, par- 
20 ticularly when reaching G new 602. This response difference is meaningful due to the listener being subjectively sensitive 
to actual linear-domain gam values during gain updates while sensitive to dB gain during steady - states of gain. That 
linear gain values tend to be of significant importance when addressing gain update issues is well understood by one 
skilled in the art and has been documented in audio literature. 

[0020] Audio muting is a special case of gain update in that the update takes place from some non-zero G old to a 
25 Gnew which is 0 Therefore, the same soft gain update may be used to implement a soft mute function. 

[0021] Figure 8 depicts the response of a constant-step soft gain update In constant-step gain update, the step- 
size, D, is the "discrete-time slope" of the line 814 joining the values of mute and full-scale or maxgain. Where mute 
- 0 and full-scale = 2 15 - 1 = 32767 in a 16-bit two's complement system a D of 8 results in 4096 steps, the closest 
integer to: 

30 

32767 32767 
N 1 ' 4096 

35 Therefore, when sweeping a gain update linearly from G old = mute to G new - full-scale, N A = 4096 gain steps are 
utilized, yielding an update time, T v of 85.33ms given a sample rate, fs, of 44.1 kHz. However, because the step-size 
is made constant, the actual number of steps for a particular gain update. N v , can differ from N-, . Figure 9 depicts the 
response of a gain update from 1/2 full scale to full-scale. If G old = full-scale/2 and G new = full-scale, only N r - N 1 /2 : 
or 2048, steps are needed. 

40 [0022] Such a constant-step approach is not as smooth as the constant-time gain update, either linear or quadratic. 
For instance, assume a 1 6-bit system with full-scale or max gain of unity (OdB). Consider two gain updates: 1 ) G okJ = 
0 = - dB (mute) and G new = maxgain/2 = 1/2 = -6dB with N,' = N 1/2 = 2048 steps; 2) G old = maxgain/2 = 1/2 = -6dB 
and G new = maxgain = 1 = OdB with N r = N 1/2 = 2048 steps. The first gain update is over nearly an 84dB range while 
the latter update is over only a 6dB range. Yet, both updates take place over the same N v and consequently the same 

45 T r where T r - N^Ts = 42.67ms. 

[0023] Figure 10 depicts in flowchart form the algorithm for controlling a linear constant-step gain update process. 
If the gam processor is busy (1002), it is currently processing a gain update and a new gain value will not be accepted. 
If the gain processor is not busy, a new value, G new , is accepted by the G host interface ( 1 004). Next if the gain processor 
has been put into bypass mode (1006), G new is immediately applied (1003). If the gain processor is not in bypass 

so mode, G new is compared to G o)d (1010), if the difference is 0, the busy flag is reset (1012) and the processor exits 
(1032) If the difference is not 0 the busy flag is set and step-counter n is set to 0 (1014) The new intermediate gain 
value G n is set to G 0 (1016) and is output as the current gain value G (1018), The value of n is incremented by D if 
G new ■ G oid > 0 or decremented by D if G new - G old < 0 (1020). When the absolute value of G new - G okJ = 0, the busy 
Nag is reset and the value of G old is set to G new (1030) and the processor exits (1032). While the absolute value of 

55 G new - G old * 0, the intermediate gain value Gn is stepped either up or down. If G new - G 0)d < 0 (1 024) then G n is stepped 
down by a number of steps, D (1026). If G new - G old > 0 1024 then G n is stepped up by a number of steps, D (1028). 
In either case, the new intermediate gain value G n is output as the current gain value G (1018). The value of n is then 
incremented by one (1020) and an exit or continuation of the gain update determination is made (1022) 



4 



EP 0 932 253 A2 



10 



15 



20 



25 



[0024] Figure 11 illustrates a gain processor which implements the above algorithm. The host interface, G Host , for 
the gain processor receives gain update signals from the host processor. In bypass mode, the new gain value G new is 
transferred through a multiplexor 1 1 04 to a register 1 1 06 holding the value G new . It is then passed through a multiplexor 
1 1 08 to a register 1110 holding the value G o)d . The value is then passed through a multiplexor 11 1 8 to a register 1 1 20. 
From register 1120, the G new value becomes the new gain output value 1124. In busy mode, the value of register 1106 
is not updated and no gain update is made. In constant-step mode, register 1114, which holds the value G n . is stepped 
in increments of D from G old to G new at adder 1116. The value of register G n is passed through a multiplexor 1118 to 
register 1120 after each step. From register 1120 the value becomes the new gain output value. A controller controls 
the bypass and busy signals of the gain processor. The controller 1124 also maintains a step counter, n and loads the 
registers 1106, 1110, and 1114. 

[0025] In constant-time gain update, the number of update steps N 1 (and consequently, the update time T^ is held 
constant regardless of the specific choices of G okJ and G new In constant-time gain updating, the step size. D, changes 
lor different pairs ol G old and G new to enable N 1 to be maintained as a constant. Figures 1 2A and 1 2B depict the nature 
of Nt in constant time updating. In Figure 12A, a gain update which sweeps from G oid - 0 to G new = maxgain takes 
steps. If = 4096. then the slope, D, 1 202 of the gain sweep = 8. Likewise, in Figure 1 2B, a gain update which sweeps 
from G old = maxgam/2 to G new = maxgain also takes N } steps However, in the latter example when N 1 = 4096, then 
the slope, D, 1202 of the gain sweep = 4. 

[0026] Generally, a constant time, variable step size gain update would require significant hardware complexity as 
A is made dependent on the gain values G oid and G new . However, the disclosed method "filters" the "instantaneous" 
gain update function to yield a first-order gain update function with a constant-time gain update. 
[0027] Figure 13 represents the linear constant-time filter output desired. The input gain change request from G okJ 
to G new merely specifies a step function. However, the disclosed gain processor will automatically translate this to a 
smoothed update, such as the linear response shown in Figure 13. 

[0028] A comb filter, having an impulse response of h(n) = 1/N-, where 0 • n*N^ - 1 and h(n) - 0 elsewhere, or having 
a z domain transfer function: 



30 



H(z)-. 



1 ) 
N i 1-z 1 



will produce a desired response curve as shown in Figure 13. When factor 1/N 1 is removed from the filter's transfer 
function, it yields the functions represented in time - domain difference equations as: 



35 



y(n)=x(n)-x(n-N 1 )+ y(n-1) 



45 



[0029] To implement these difference equations, the term x(n - N,). which, for ~ 4096 as an example, becomes 
x(n - 4096), would require a storage of 4096 input signal values to be obtained. However, by observing the fact that 
input signal x(n) takes on only two values, G oid and G new , and has a step-function profile as depicted in Figure 13, 
simple external control circuitry can be designed to generate the state variable x(n - Nj). 

[0030] Figure 14 schematically depicts the implementation of the above difference equation. The initialization and 
values, when bypass = 0, during update of the registers x n , x n . N1 , and y n _j are as follows: 



50 



Counter 


x n 


X n-N1 


y n -i 


9n 


n = 0 


G old 


G old 


Gold* 


G old 


0 < n <N 1 - 1 


Gnew 


G old 


y n -i 




n = Nt 


new 


°new 


G new* N, 


Gnew 



[0031] The host interface, G Host , receives a gam update command from its host. Gain controller maintains the busy 
ss and bypass indicators for the gain processor. The gain controller also maintains an internal counter to track the sweep 
update and to maintain the busy mode. The gain controller also loads values into registers 1410, 1412, and 1418, 
which hold the values for x n x n . N1 and y n .,. respectively Registers 1 410 and 141 2 initially hold the value G old Once 
the counter maintained by gain controller increments past 0, the gain controller loads register 1410 with value G new . 
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The value of register 1 41 2 does not change during the gain update process. Adder 1 41 4-subtracts the value of register 
1412 from the value of register 1410. Adder 1416 adds the value of register 1418, which is initially loaded with value 
G oid* N i' to tne output of adder 1414. Also, once the counter maintained by gain controller increments past 0. the gain 
controller 1404 loads register 1418 with value y n . v When the counter ts equal to the number of steps, N 1( register 1418 
holds the value G new *N 1 . Value y n takes its value from the result of adder 1416. Shift 1422 is a right shifter which 
divides value y n by N } to produce the output G n . 

[0032] The counter maintained by the gain controller may take on values of 0 to N, - 1 only. The case of n = N, may 
be interpreted as the counter rolling over to 0 with Gold being replaced by G new . The gain controller, during 0 < n < N 
- 1 would be in the busy mode and therefore would not accept any new gain updates from the host 
[0033] For = 4096 and an output gain, g(n), of 16-bit gain, the addition operation, x(n) - x(n - 1 ), requires an adder 
141 4 of 17-bits as registers 1410 and 1412 would be 16 bits each. The feedback addition operation requires an adder 
1416 of 30 bits as register 1418 is 30 bits in size. 

[0034] The algorithm for linear, first-order, constant-time updates described above is extended to achieve a quadratic, 
second-order, update. In this case, a second-order filter response as depicted in Figure 1 is desired. The input gain 
change request from G o)d to G new merely specifies a step function. However, the disclosed gain processor will auto- 
matically generate a smoothed update, such as the quadratic response depicted in Figure 1 . The response may be 
obtained using a second-order comb filter characterized by M - N 2 /2. where M represents the midpoint of the number 
of steps in the update. 

[0035] The polynomial fit is carried out on two time ranges: 0 < n < N and N < n < 2N. The midpoint of the curve, N, 
ts defined by {G new + G o)d )/2- For values on the curve < N , g(n) = c^n 2 + a, where g(0) = G okj and g(N) = G new . Therefore, 
= ( G new * G old )/2 and a, = G okJ . For values on the curve > N and < 2N, g(n) = P 0 (n - 2N) 2 + p, where g(N) = G new - 
G old and g(2N) = G new . Therefore, p 0 = -(G new - G old )/2 or - p 0 and ^ = G new . 

[0036] The time-domain transfer function of the comb filter which yields the response depicted in Figure 1 is: 

25 

y(n)- x(n)-2x(n-M)+x(n-2M)+2y(n- 1 )-y(n-2) 



30 




A sample implementation of the above time-domain transfer function as a comb filter is described in the following 
matiab code: 

35 

comb2 .m 

40 % implementation of the 2nd-order comb filter 
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% integrator is implemented as a second-order section 
% 

% Prashant P. Gandhi, August 1997 
% Algorithm Implementation Code 
% 

% H(z) * M * M = [1-2* z A (-M) + z^(-2 * M)]/[l - 2 *2^(-l) + 
z A (-2)J 

% y(n) = x(n) - 2 * x(n-M) + x(n - 2 * M) + 2 * y(n - 1) - y(n 

- 2) 
% 

clear all; 
globdef ; 
regdef ; 



N = 11;' M = 4; GO = 0; Gl = 2~(gainbits - 1) - 1; 

bb = [1 zeros (1,M - 1) - 2 zeros (1,M - 1) 1]; aa = [1 - 2 1] ; 

x = [Gl Gl*ones(l,N - 1)] 



xy = filter (bb, l,x) ; 
y = f ilter (1, aa, xy) ; 
ym = y/ (M * M) ; 
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40 
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xn = 0; xnM = 1; xn2M = 2; 

ynlL = 3; ynlH = 4; yn2L = 5; yn2H = 6 ; yn = 7; x2tmp = 8, 
yLtmp = 8; 



for n = 1:N; 
if rem (n, 25) = =0, fprintf('%d \n); end; 
if n == 1 
yoldl(n) = 0; yold2(n) = 0; 

elseif n == 2 

yoldl(n) = y2 (n - 1); yold2 (n) = 0 
else 

yoldl(n) - y2(n - 1); yold2(n) - y2(n - 2) 
end 

if n< = M 

55 xoldM(n) = GO; xold2M(n) = GO 



7 



elseif n> = M+l & n< 

xoldM(n) = Gl; 
else 

xoldM(n) = Gl; 
end; 
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= 2 * M 

xold2M(n) = GO 

xold2M(n) = Gl; 



xl(n) = x(n) -2 * xoldM(n) ; 

x2(n) = xl(n) + xold2M(n) ; 

yl(n) = x2(n) + 2 * yoldl(n); 

y2(n) = yl{n) - yold2 (n) ; 

y3(n) - y2(n)/(M * M) ; 



% Microcode DSP Instructions 
% 

E = x{n) ; 

0; nopalu([]); st(e, xn) ; 

xoldM(n); 
0; st(e,xnM); 
E = xold2M(n) ; 
0; 



i - 



xx (n) = E; 



i = 

i = 

i = 

i = 



1; 
2; 

3; 

4; 



unsigned 
i = 5; 
i = 6; 



carryl 
i = 7; 

i = 9; 

carry2 



st (e, xn2M) 
ld(xnM,bc) 
add(c,b, 0,b) ; ld(xn,a) ; 



xxnM(n) = B; 
xxn(n) = A; 
%B = 2 * x<n - M) ; 
sub(a,b, 0,b) ; ld{xn2M,a); xxn2M(n) = A; 

xxl(n) = B; %B = xl(n) 

add(a,b, 0,a) ; ld(ynlL,bc) ; xx2 (n) = A; 

yoldlL(n) = B; %A = x2(n), xx2(n) 



%save x2 



is 



st (a, x2tmp) ; 
addgO (c,b, 0, a) ; ld(ynlH,bc) ; yoldH(n) = B; 

zyoldlL(n) = A; carryl(n) = c2out_reg; 

%A = 2 * ynlL, gen 

addgl (c,b, 0, c) ; ld(x2tmp,b) ; zyoldH(n) = C; 

%C = 2 * ynlH + carryl 
addgO (a,b, 0, a) ; ylL(n) = A; carry2 (n) = c20out_reg; ' 

%A = 2 * ynlL + x2 , gen 
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i = 11; 



12; 

13, 
14; 
15; 
16; 
17; 
18; 
19; 



addgl<c,bmsb,0,c) ; ld(yn2L,b); ylH(n) =C; 

yold2L(n) = B; %C = 2 * ynlH + carryl + carry2 
subg0(a,b, 0,a) ;ld(yn2H,b) ; yL(n) =A; 

yold2H(n) = B; carry3 (n) = c20out_reg; 

%A = ylL - yn2L , gen carry3 
subgl(c,b, 0,c) ; st (a,yLtmp) ; yH(n) = C; 

%C = ylH - yn2H + carry3 
ld(ynlL,a) ; 
st (a,yn2L) ; 
ld(ynlH,a); 
st (a,yn2H) ; 
ld(yLtmp, a) ; 
st (a,ynlL) ; 
st(c,ynlH); 



% Check for error between algorithm implementation and 
microcode 

30 % 

err = sum(abs{yL, + yH * M * M - y2 ) ) ; 
if err- = 0 

35 fprintf('\n ERROR in Gain Computation^ ' ) 

else 

fprintf('\n Gain Computation is OK\n') 
end; 

40 

[0037] For M = 1 024, x(n - M) and x(n - 2M) would require storage of 2048 input signal values. However, the algorithm 
shown in the table below allows an external controller that generates state variables to eliminate the 4 kilobyte word 
storage requirement for a 16-bit gain value. 

45 



Cond. and Count 




X n-M 


X n-2M 


y n -i 


yn-2 


9n 


G old = G new n = 0 


G old 


G old 


G old 


Goid'M 2 




G old 


G old= G new n = 0 


G new 


G old 


G old 


Goid'M 2 


Goid'M 2 


G ne w/M2 


1 <n<M- 1 


G new 


G old 


G old 




y n . 2 


y n /M2 


M < n < 2M - 1 


°now 


G new 


G old 




y n - 2 


y n /M2 


n >2M 


G new 


G new 


G new 


GnewM 2 


G„ewM 2 


G new 



For n > 2M, the counter is reset to 0 with G dd = G new , and G new = G host . 

[0038] Figure 1 5 depicts a sample ol a gam processor which implements the above algorithm table, where M = 1024 
or N 2 = 2048. Register 1502, which holds the value x n , is initially loaded with value G new . Registers 1504 and 1506, 
whtch hold the values x n . M and x n . 2M , respectively, are initially loaded with the value G old . The value of y n-1 and y n . 2 at 
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the first step, n = 0, is G otd * M 2 . The output, g(n) is G okJ . Then, while the counter, n, is less than or equal to 1/2 of the 
steps in the gain sweep, 2M, registers 1502, 1504, and 1506 remain unchanged. Multiplier 1508 multiplies the value 
of register 1504 by two. Adder 1510 subtracts the product of multiplier 1508 from the value of register 1502. Adder 
1512 adds the result of adder 1510 to the value of register 1506 The value of y^, is dictated by the feedback loop of 

s y(n) and unit delay 1514. The value of y n0 is multiplied by two at multiplier 1518. The value of y n _ 2 is generated by the 
unit delay 1520 tn the feedback loop of y(n). Adder 1524 subtracts the value of y n . 2 from the product of multiplier 1518. 
Adder 1 526 adds the result of adder 1 51 2 to the result of adder 1 524 The result of adder 1 526 is y(n) This result is 
input into Shift 1 530 a right shifter which divides the value y n by M 2 to produce the output gain signal, g(n). The counter, 
n, is incremented after each series of operations. While the counter, n, is greater than or equal -to 1/2 of the steps in 

10 the gain sweep, 2M t register 1504 is set to G new . The above steps remain unchanged. Once the counter, n, reaches 
2M. registers 1502, 1504, and 1606 and g(n) are valued at G new . 

[0039] A controller for the gain processor depicted in Figure 15 preferably implements the following features: apply 
a gain of 0 (mute) on power up, maintain a busy state during gam updates to prevent gain from being altered; while in 
a steady-state, the applied gain would be forced to equal G ofd allowing for recovery from "rare" gain errors which could 
*5 be caused by noise on a line, faulty power supplies, or errors in the physical device; maintain a bypass state such that 
G oid r G new = Ghost with the applied gain being G old . 

[0040] Figure 16 depicts a microcode cycle-based DSP implementation of the table algorithm with a single adder 
The filter may be efficiently implemented, for an example, with g(n) represented in 20 bits and M = 1024, that is, N 2 = 
2048 or T 2 = 46 44ms at fs = 44 1kHz The RAM words of the DSP are assigned, in this embodiment, as follows: 

20 [0041] G 0 | d , G new = x n , x n . M and x n . 2M : lower 20 of the 32-bits with gain values and upper 1 2 bits set to 0. 

[0042] The values y n , y n _-,, and y n . 2 are assigned two RAM words each: least significant word, LSW, and most sig- 
nificant word, MSW. The lower 20 bits of the MSW of y n , y n0 and y n . 2 are set to G o)d while the remaining upper 1 2 bits 
as well as all 32 bits of the LSW of y n , y^ and y n . 2 are set to 0, during initialization. 
[0043] The MSW of y n = g n with the lower 20 bits being the valid g n . 

25 [0044] A multiplexor 1602, controlled by a gain controller, is used as a data selector to provide input to the RAM 
words described above. Registers 1604, 1606, and 1608 are 32 bit registers which are loaded with values, x^ x n . M , 
X n-2M> Yn> V n-1> and y n . 2 from RAM as needed to perform the appropriate step of the filter. A multiplexor 1610 is used 
as a data selector, determining which value to load into the adder 1612 and 1616. The two feed-forward additions (x 
(n) + (-2x(n - M)) + x(n - 2M) are carried out as normal additions Each feedback addition 1 524 and 1 526, on the other 

30 hand, is carried out in two clock cycles: one cycle to add LSWs and another cycle to add the MSWs. Special addition/ 
subtraction instructions are used for LSW and MSW operations. When the LSW add is performed, the Cout of the 
20-bit adder 1616 is stored into flip-flop C20. The value of flip-flop C20 is used as the Cin to the 20-bit adder 1616 
during the MSW addition operation. The value, 2+y^ is implemented by adding two LSWs of y n . 1( storing the Cout of 
adder 1616 in C20, passing it as the Cin of the 12-bit adder 1612 at the next clock and adding the two MSWs of y n . v 

35 Shifter 1628 carries out the multiply by two operations of multipliers 1508 and 1518 and the shift function of 130. A 
register 1626 holds the value g(n) for output. 

[0045] In the presently preferred embodiment, this DSP is implemented with a datapath size which is less than 40 
bits. This 32 bit implementation allows for a 25-33% processor area reduction. However, twice as many addition op- 
erations and twice as many RAM storage words are required for feedback processing. If DSP area cost is not a concern. 

•*o the above implementation may be altered to accommodate a 40-bit or larger datapath. 

[0046] Figure 2 depicts a block diagram of the signal processing path of a sample audio system which advantageously 
utilizes the soft gain update. An equalizer, audio tone control, and volume control are each part of a digital signal 
processing hardware which also includes a DSP control processor, control circuitry 218 and memory. The equalizer 
receives digital audio data. The data is processed by the DSP hardware. Digital to analog circuitry converts the signal 

^5 for output as analog audio. When a gain update command is received by the DSP control processor from user input, 
the DSP control processor commands the volume control block, via the DSP control circuitry, to update the gain ac- 
cordingly. The volume control block utilizes.the soft gain update. 

[0047] In a sample implementation of linear constant-time soft gam updates, the best mode presently contemplated 
for sweep time is in the neighborhood of 90ms. In a sample implementation of quadratic constant-time soft gain updates, 

50 the best mode presently contemplated for sweep time is in the neighborhood of 45ms. Linear sweep times of between 
70ms and 11 0ms and quadratic sweep times of between 35ms and 55ms are preferred because they provide the 
listener a pleasant gain update response. Linear sweep times of less than 70ms or more than 110ms and quadratic 
sweep times of less than 35ms or more than 55ms are less preferred due to a greater amount of audible artifacts 
occurring when such sweep times are used. 

55 [0048] As will be recognized by those skilled in the art, the innovative concepts described in the present application 
can be modified and varied over a tremendous range of applications, and accordingly the scope of patented subject 
matter is not limited by any of the specific exemplary teachings given. 

[0049] In the preferred embodiment, a 16-bit audio application is contemplated. However, audio applications of var- 
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ying bit resolutions can take advantage of the present disclosure. 

[0050] The number of steps between mute and maxgain in the preferred embodiment is 4096. (3276S levels with a 
slope, D, of 8). However, the slope can be varied, by either powers of two or some other number, to produce a steeper 
or flatter response. 

[0051] Like the slope, the sample rate, fs, can differ among audio applications. As the sample rate changes, the 
sweep time for constant-step soft gam update wilt vary accordingly 

[0052] The soft gain update described above is not limited only to first and second order functions. It can be used 
with higher ordered functions as well, each function having a greater continuity at G new and G old . Therefore, each 
higher ordered functions will produce fewer audible artifacts at gam update. However, implementation of such higher 
ordered function will demand a greater amount of DSP resources Preferable sweep times for higher ordered function 
are in the neighborhood of one-half of the preferred sweep time for the function of the next lowest order. That is, a 
preferred third order soft gain update sweep time will be in the neighborhood of 22.5ms (1/2 of 45ms, the preferred 
sweep time for a quadratic update). 

[0053] It should also be noted that the three main embodiments above are not exhaustive, and other algorithms can 
be used for smooth automatic updating. 



Claims 



20 1. An audio amplification system, comprising: 
audio amplifier circuitry; and 

a controller for controlling the gain of said amplifier circuitry, and for changing the gain of said amplifier circuitry 
according to a predetermined romping function in response to a gain-change command. 

25 

2. The system as claimed in Claim 1 , wherein the said controller is arranged for changing the gain of the amplifier 
circuitry in a plurality of steps. 

3. The system as claimed in Claim 2, wherein at least one of said plurality steps occurs after the end of said gain- 
30 change command. 

4. The system as claimed in any preceding claim, wherein the controller is arranged for changing the gain of the 
amplifier smoothly. 

5. The system of Claim 4 wherein the controller is arranged for changing the gain of the amplifier over a predetermined 
and extended period. 

6. The system as claimed in any preceding claim, wherein the controller is arranged for changing the gain of the 
amplifier according to a predetermined smoothing function of at least the second order. 

40 

7. The system as claimed in any preceding claim, wherein said controller is integrated with said amplifier circuitry 

8. The system as claimed in any preceding claim, wherein said controller comprises car-radio control functions. 

is 9. The system as claimed in any preceding claim, wherein said controller comprises a DSP core that is integrated 
with custom peripherals for I/O functions. 

10. The system as claimed in any preceding claim, wherein said controller is arranged for changing the gain in less 
than 110 milliseconds. 

so 

11. The system as claimed in any preceding claim, wherein said controller is arranged for changing the gain in less 
than 55 milliseconds 

12. The system as claimed in any preceding claim, wherein said controller is arranged such that it does not respond 
55 to successive gain-change commands while it is changing the gain of said amplifier circuitry. 

13. The system as claimed in any preceding claim, wherein said controller is arranged for changing the gain of said 
amplifier circuitry automatically. 
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14. A method of gain update of a signal comprising: 

receiving a gain-change command; 

automatically changing the gain of said signal in a plurality of steps according to a predetermined ramping 
5 function. 

1 5. The method as claimed in Claim 1 4 further comprising; 
performing the step of changing the gain in less than 110 milliseconds. 

io 1 6. The method as claimed in Claim 1 4 or 15 further comprising: 

preventing performance of said step of receiving the gain change command while the step of automatically chang- 
ing the gain is performed. 

is 
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